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Abstract: [ Aim] Reticulitermes flaviceps is one of the representative wood-feeding termite in our country, 


which has efficient cellulolitic system. This study aims to clone and characterize cellulase gene from R. 


flaviceps, so as to provide the theoretical basis for exploring cellulose resource of termites. [Methods] A B- 


glycosidase gene was cloned from R. flaviceps by using reverse-transcription polymerase chain reaction ( RT- 
PCR) , degenerate PCR and rapid amplification of cDNA end (RACE) methods, and its sequence was analyzed 


by bioinformatics tools. The corresponding protein was prokaryoticly expressed in Escherichia coli and the 


cellulase activity was measured after purification. 【 Results] The full-length cDNA sequence of B-glucosidase 
gene GZRfbg| (GenBank accession no. ; KU886065 ) was cloned from R. flaviceps, with a length of 1 691 bp 


and an open reading frame of 1 488 bp, encoding 448 amino acid residues with the typical structure of glycosyl 


hydrolase family (GHF) 1 domain. The calculated protein molecular weight is 56.45 kD. Multiple alignment 
analysis revealed that the deduced amino acid sequence of GZRfBG1 shares 95% identity with the salivary gland 
B-glucosidase of Reticulitermes flavipes. GZR{BG1 was heterologously expressed in Æ. coli Rosetta 2 (DE3). 
The purified tagged GZRfBG1 exhibited increased specific activity to salicine, carboxymethyl cellulose sodium 


(CMC) , and cellobiose. And the cellulolytic activity was also shown in a CMC containing agarose plate, in 


which clear digestion haloes were formed. [ Conclusion] These results suggested that the salivary gland B- 


glucosidase from sibling species R. flaviceps and R. flavipes are very conservative in sequence and hydrolysis 


characteristics. This study of B-glucosidase gene can enrich the theoretical basis for the development of termite 


cellulases and phylogenetic study. 
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1 INTRODUCTION 


Lower wood-feeding termites, with profoundly 


efficient abilities to digest cellulose, play an 


important role in the recycling of wood 
lignocelluloses (Feng et al., 2015; Peterson and 
Scharf, 2016 ). 


cellulose 


These termites possess a dual 


degradation mechanism. Endogenous 
cellulases and symbiotic microorganism ( bacteria 
and protozoa) cellulases cooperatively contribute to 
cellulose digestion. This system has three main types 
of cellulases, i. e. , endoglucanases (EC 3.2.1.4), 
cellobiohydrolase ( EC 3. 2. 1. 91), and B- 
glucosidases (EC 3.2.1.21) (Scharf et al., 2011a, 
20116; 2014 ). 


hydrolyze glucose dimers and cello-oligosaccharides 


Brune, B-Glucosidases can 


to glucose, and this is essential to complete cellulose 
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digestion (Tarayre et al., 2015). Lower termites can 
secrete endogenous Q-glucosidases in their salivary 
glands and midguts ( Watanabe, 1998; Lo and 
Watanabe, 2011). The salivary gland-specific B- 
glucosidase has been analyzed via heterologous 
species , 
including Neotermes koshunensis ( Shiraki) (Ni et 
al., 2007; Uchima et al., 2011), Nasutitermes 
takasagoensis Shiraki (Tokuda et al., 2002, 2009; 
Uchima et al., 2013 ), Coptotermes formosanus 
(Shiraki ) (Zhang et al., 2010; Zhang et al., 
2012a ) and Reticulitermes flavipes ( Koller ) 
(Hirayama et al., 2010; Scharf et al., 2010). 


Recombinant cellulolytic enzymes are useful for 


expression systems in several termite 


studying cellulose metabolism in termites and have 
the potential to be engineered as biocatalysts for 
cellulose-based biofuel production. 
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The wood-feeding termite Reticulitermes flaviceps 
(Oshima) is widely distributed in southern China. 
The salivary gland B-glucosidase of R. flaviceps was 
reported firstly in 2008. Its nucleotide or amino acid 
sequence remains unknown, but its estimated 
molecular mass is 93.6 kD. In this study, we cloned 
an endogenous B-glycosidase gene named GZRfbg1 
( glycoside hydrolase family 1 gene from the salivary 
gland of R. flaviceps). Its encoded protein has the 
molecular mass of ~ 54. 6 kD. Bioinformatics 
analysis and heterologous expression of GZRfbg1 
were then performed. We expect that these results 
understanding of the 


will enrich our termite 


lignocellulase system. 


2 MATERIALS AND METHODS 


2.1 Termites 

Workers of R. flaviceps were collected from a 
laboratory maintained single colony and washed with 
double distilled water before use. 
2.2 cDNA cloning of B-glucosidase gene from 
R. flaviceps 

Total RNA and cDNA of the whole termite body 
were obtained using a Total Isolation System Kit and 
the PrimeScript™ II 1st Strand cDNA Synthesis Kit 
(TaKaRa, Japan) , respectively, using 1.2 ug RNA 
as a template for cDNA synthesis. To obtain cDNA 
fragments of the target gene, PCR was carried out 
using degenerate oligonucleotide primers that were 
synthesized based on the other identified partial 
amino acid sequences of B-glucosidases, and the 
primers used for degenerate PCR were designed by 


the iCODEHOP program (http ://blocks. fherc. org/ 


codehop. html ) and shown in Table 1. The 
degenerate PCR reaction system of 50 uL contained 
each primer (2 pmol/L) 2 uL, 25 uL of Taq DNA 
polymerase mixture (TaKaRa, Japan). The thermal 
cycling conditions were 94°C for 2 min, followed by 
30 cycles of 94 for 30 s, 50 for 30 s, 72°C for 
30 s, and then 72°C for 10 min. Based on the 
sequences of degenerate PCR products, primers 
shown in Table 2 were designed for 3'-RACE and 5’- 
RACE. One round 3'-RACE, three rounds 5'-RACE 
and nested-PCR were 
manufacturer’ s instructions (Invitrogen, GeneRacer™ 
Kit, USA ). The 


sequenced and assembled. The complete cDNA was 


performed according to 


amplified fragments were 
amplified and cloned into a plasmid vector pMD-20T 


(TaKaRa, Japan). 
selected for sequence analysis and elimination of the 


Five clones were randomly 


PCR mutations. The final determined sequence was 
examined using the BLAST program provided by 
NCBI (http ;//blast. ncbi. nlm. nih. gov/Blast. cgi). 
The cleavage site of the signal peptide was predicted 
by the SignalP program (http://www. cbs. dtu. dk/ 
services/SignalP ). 
predicted using a combination of neural networks 
(NN) and hidden Markov models (HMM) (http :// 
www. cbs. dtu. dk/services/NetNES/). The NetPhos 


2.0 server neural network program was used to predict 


Nuclear export signals were 


serine, threonine and tyrosine phosphorylation sites 
( http://www. cbs. dtu. dk/services/NetPhos/ ) . 
Prediction of transmembrane helices in proteins was 
accomplished using the transmembrane helices 
program on the TMHMM Server v. 2.0 (http:// 
www. cbs. dtu. dk/services/TMHMM/ ) . 


Table 1 Degenerate primers used to amplify the B-glucosidase gene from Reticulitermes flaviceps 


GenBank accession numbers of the corresponding 





Primers Primet sequences (9 = 3) homologous protein sequences 
Fl CCAAAGGAGTTAATATHTGGGAYAC 
F2 TACCCCGGGGATCGGNGAYTAYYT 
F3 ACTTCTTTGGCATGAAYTWYTAYAC 
R1 TANATRGTRACCCTGGACGGCGT gi269965728, gi269965726, gi269965730 ; 
R2 ATRGTRACCCTGGACGGCGTCTTTG 21269965724, gi269965718, gi269965720 ; 
R3 CCNYTRATRGGAAGGCATCACTTTC gi291293226 ，gi303324839 , gi303324840 ; 
R4 CTYCCNAWRTGATGGAGAGCCGAT gi300432455 ,gi283484500 
R5 TTRACCTADYKTTTGCTCATACCCT 
R6 ATRATRTCNCTCGTCTACTTCCTCT 
R7 CTRTTRAARCTCACCGACGCTCCC 


Degenerate nucleotides are as follows; R= A/G, Y =C/T, M=A/C, K=G/T, S=C/G, W=A/T, H=A/C/T, B=C/G/T, V=A/C/G, D=A/G/ 


T, N=A/C/G/T. 
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Table 2 Primers used in RACE 


Primers 


Primer sequences (5’ —3’) 


Use of primers 





R. BIB7-3'Race-sp1 


R. B1B7-3'Race-sp2 


CGGAATTCCAGTGTTTTGAAAGCGAGGGTAGAC 


CGGAATTCGCCTGGATTGAGCCGGTGACTAA 


3'-RACE 


AP GGCCACGCGTCGACTAGTACT TTTTTTTTTTTTTTTT 


AUAP GGCCACGCGTCGACTAGTAC 





R. BIB7-5'Race-sp1 
R. BI B7-5'Race-sp2 
5’Race abridged* 


5'Race abridged anchor 


CGGAATCCCCATGGTACAACTCTG 


TGAATTTTGGTAGACGGGATGTGG 


5’-RACE Round 1 


GGCCACGCGTCGACTAGTACGGGHGGGIGGGHG 


GGCCACGCGTCGACTAGTAC 





R. B1B7-5'Race-sp4 


B1B7-5'Race-sp5 


CTCGCTTTCAAAACACTGGGGTAA 


5’-RACE Round 2 


AACCTCTCCTTTCTGATCCGCTTTGA 





B1B7-5'Race-sp6 


The “II” in 5’-RACE abridged primer represents poly “G”. 


2. 3 Expression of recombinant protein in 
Escherichia coli 

The sequence of -glucosidase gene, without 
the nucleotide sequences encoding the putative N- 
terminal signal peptide, was amplified by PCR and 
analyzed by sequencing. The PCR primer sequences 
were as follows: the forward primer was flanked at 
the 5’ end with an Ncol restriction sequence (5’- 
CCATGGGACGTCGATAACGAAACCCT-3' ) and 
the reverse primer flanked at the 5’ end with an Xhol 
restriction sequence omitting the stop codon (5’- 
CTCGAGGTCTAGGAAGCGTTCTGGAAT-3') for 
His-tag fusion expression. The underlined was used 
to indicate the recognition site of restriction enzyme. 
The PCR reaction system of 50 uL contained a final 
concentration of 0. 1 pmol/L of each primer, DNA 
templates at a final concentration of 100 ng, and 25 
uL of Taq DNA polymerase Mixture ( TaKaRa, 
Japan). The thermal cycling conditions were 94°C 
for 2 min, followed by 30 cycles of 94°C for 30 s, 
55°C for 30s, 72°C for 1 min, and then 72°C for 10 
min. The cellulase gene overlapped with Nco I and 
Xho I cohesive ends was cloned into pET32a vector 
(Novagen, USA) at corresponding restriction sites. 
Recombinant pET32a plasmids were propagated in 
Escherichia coli DHS a and five clones were randomly 
selected for sequencing verification. Competent cells 
of E. coli BL21 and Rosetta 2 (DE3) (Novagen, 
USA ) were individually with the 
following 2 plasmid constructs; ( 1) recombinant 


transformed 


pET32a plasmid carrying the cellulase gene without 
the stop codon at the 3’ end, thereby allowing the 
expression of recombinants with 6 x His residues 
downstream of the gene (referred to as tGZRfbg1 
and (2) 


plasmid with no cleavage or insertions (referred to as 


hereafter ) , original control pET32a 


GTAACCATCGGCTCTATTTCATTGTC 


5’-RACE Round 3 


pET32a-CK hereafter ) . 
randomly picked for growth in Luria-Bertani ( LB) 


Resulting colonies were 


medium with ampicillin (100 pg/mL) in shaker 
flask. E. coli cells were incubated at 37°C at 220 1/ 
min and induced with 0.5 mmol/L isopropyl B-D-1- 
thiogalactopyranoside ( IPTG ) ( Dingguo 
Changsheng, China) for an additional 3 h when the 
optimal density reached 0. 4 - 0. 5 at 600 nm. 
Induced E. coli cells were divided into two aliquots 
(A and B) and harvested by centrifugation; (A) 
the cell pellet was suspended in 50 mmol/L Tris-HCl 
buffer (pH 8.0) containing 5 mmol/L EDTA and 1 
mmol/L PMSF ( Dingguo Changsheng, China). The 
cells were lysed by sonication in ice-water and 
centrifuged at 8 000 x g at 4°C for 20 min. The 
supernatant and the insoluble matter fractions were 
boiled for 5 min after adding 2 x SDS-PAGE (sodium 
dodecyl sulfate polyacrylamide gel electrophoresis ) 
loading buffer to detect soluble protein and proteins 
expressed as inclusion bodies (IB), respectively; 
(B) the cell pellet was suspended in 2 x SDS-PAGE 
loading buffer and boiled for 5 min then used for 
analysis of total cell protein (TCP). SDS-PAGE and 
Western blot were used to identify the expression of 
tGZRfbg| and carried out referring to Molecular 
Cloning II (Sambrook, 2001). The protein bands 
difluoride 


(PVDF) membrane and successively incubated with 


were blotted onto a polyvinylidene 
primary antibody ( mouse-anti-His-tag ) and HRP- 
conjugated anti-mouse IgG (MBL, Japan). The 
immunoreactive bands were visualized with a DAB 
(3, 3’ N-diaminobenzidine — tertrahydrochloride ) 
Western Blotting Detection Kit (CWBIO, China). 
2.4 Isolation and purification of recombinant 
protein 

For large-scale isolation and purification of 
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recombinant protein, recombinant E. coli Rosetta 2 
(DE3 ) cells were grown in 100 mL LB with 
antibiotics as described in section 2. 3 and shake- 
incubated at 37°C for 2 -3 h at 220 r/min and 
induced with 0.1 mmol/L IPTG for an additional 18 
h at 20°C when the optimal density reached 0. 4 at 
600 nm. The cells were harvested by centrifugation 
and lysed by the sonication method described in 
section 2. 3. The supernatant was recovered by 
centrifugation at 5 000 x g at 4°C for 20 min and the 
cell debris was removed by passage through a 
membrane filter (中 =0.45 um, Millipore, USA). 
The purification of soluble tGZRf{BG1 was performed 
according to the manufacturer’ s instructions using a 
6 x His-tagged ( Ni-Agarose ) Soluble Protein 
( CWBIO, China ). 
purification, the purified protein was stored in 100 
mmol/L sodium acetate (SAB), pH 5.6, with 10% 
glycerol, using an Ultra-15K centrifugal filter unit 
(Millipore, Cork, Ireland). Purified tGZRfBG1 was 
preserved in SAB (pH 5.6) with 30% glycerol. 
Protein concentration was estimated with the 
Bradford Reagent (CWBIO, China) using bovine 
serum albumin (BSA) as the standard protein. SDS- 
PAGE and Western blotting were used to identify the 
purified tGZRfBG1 and performed as previously 


Purification Kit Following 


described in section 2.3. 
2.5 Activity assays of recombinant cellulase 
The pET32a-CK, expressed and purified under 
the same conditions as tGZRfBG1 in section 2. 4, 
was used as the negative control protein. The 
cellulolytic activity of purified soluble tGZRfBG1 was 
evaluated by the following assays: (1) Agarose plate 
assay: the agarose plate was made with 1% agarose 
and 0.2% carboxymethyl cellulose sodium (CMC) 
The purified tGZRfBGI and 
pET32a-CK solutions were spotted on the surface of 


in distilled water. 
the plate in equal amounts. The plate was then 
incubated at 37°C for 18 h. The enzyme activity was 
visualized using the Congo Red Stain method ( Zhang 
et al., 2009). (2) Substrate specificity assay: the 
enzyme was incubated with 1% CMC, 1% salicine, 
and 1% cellobiose, respectively. Enzyme activity 
was assayed by measuring the release of glucose from 
this substrate using the dinitrosalicylic acid (DNS) 
method (Zeng, 2016). Each assay was performed in 
triplicate. One unit (U) of specific enzyme activity 
was defined as the amount ( pmol) of reducing 
sugars (glucose equivalent) produced by 1 mg of 
protein per min. An equal amount of pET32a-CK 
solution was used as the negative control. 


3 RESULTS 


3.1 Sequence features of B-glucosidase gene of 


R. flavicepes 

The full-length cDNA of B-glucosidase gene of 
R. flavicepes obtained in this study was 1 691 bp and 
designated as GZRfbg1 (GenBank accession no. : 
KU886065). This cDNA contains an open reading 
frame (ORF) of 1 488 bp. The translation initiation 
codon ATG in cDNA nucleotide position is from 45 to 
47 bp and the termination codon TAA is from 1 530 
to 1 532 bp. In the 3’-terminal region, there is a 
putative poly-adenylation signal sequence AATAAA 
and a poly (A + ) tail, suggesting that the GZRfbg1 
was derived from an eukaryotic organism (Fig. 1). 
The ORF encodes 448 amino acids with the 
calculated molecular weight of 56. 45 kD. This 
multiple alignment results showed that GZRfBG1 has 
the typical structure of glycosyl hydrolase family 
(GHF) 1 domain and shares 95% , 87% , 86% , 
83% and 81% amino acid sequence identity with the 
salivary gland B-glucosidase sequences of R. flavipes 
( GenBank accession no.: ADK12988.1 ), 
Coptotermes gestroi ( GenBank accession no. : 
AGS32242.1), C. formosanus (GenBank accession 
no. : ADB23476. 1 ) Macrotermes barneyi Light 
AFD33364. 1 ) and 
Nasutitermes takasagoensis ( GenBank accession 
no. : BAI50023. 1) , respectively (Fig. 2). A signal 
peptide cleavage site between the 17th and 18th 
amino acid 


( GenBank accession no. : 


Generic 
phosphorylation sites analysis results showed that 
GZRfBG1 has 12 serine, 4 threonine and 9 tyrosine 
phosphorylation sites. A nuclear export signal 
(NES) was detected at the 430th amino acid 
residue. Transmembrane helices analysis indicated 
that GZRf{BG1 is probably located outside the 
membrane. All the forecasting mentioned above 
suggested that GZRfBG1 has 


modification and can be secreted from cells into the 


residue was found. 


post-translational 


surrounding tissues. 
3. 2 Expression and purification of the 
recombinant GZRfBG1 

Over-expressions of tGZRfbgl and pET-32a 
vector ( pET32a-CK ) were detected by IPTG 
induction of LB culture, and analyzed by SDS-PAGE 
and Western blot. The expression of tGZRfBG1 was 
not detected in E. coli BL21 cells. Most tGZRfBG1 
was expressed as inclusion bodies in Æ. coli Rosetta 
2 (DE3) cells. The expected expressed fusion 
protein encoded by pET32a-CK alone would be 
approximately 26 kD. The tGZRfBG1 was expressed 
with Trx. tag, S. tag and His-tag. In comparison with 
of pET32a-CK ( negative 
control ), a distinctive band with calculated 
molecular weight (MW) around 74 kD was found in 


the protein profile 
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Poly-adenylaton signal sequence 


AATAAA 


polyA-3’ 


ORF: 1 488 bp 


Tnitiation codon ATG Termination codon TAA 
45-47 bp 1 530-1 532 bp 


cDNA: 1 691 bp 











Fig. 1 Schematic structure of GZRfbg1 cDNA of Reticulitermes flaviceps 





Reticulitermes flaviceps 94 
Reticulitermes flavipes 94 
Coptotermes gestroi 94 
Coptotermes formosanus 94 
Nasutitermes takasagoensis 91 
Macrotermes barneyi 94 
Reticulitermes flaviceps 188 
Reticulitermes flavipes 188 
Coptotermes gestroi 188 
Coptotermes formosanus 188 
Nasutitermes takasagoensis 185 
Macrotermes barneyi 188 
Reticulitermes flaviceps 282 
Reticulitermes flavipes 282 
Coptotermes gestroi 282 
Coptotermes formosanus 282 
Nasutitermes takasagoensis 277 
Macrotermes barneyi 260 
Reticulitermes flaviceps 376 
Reticulitermes flavipes 376 
Coptotermes gestroi 376 
Coptotermes formosanus 376 
Nasutitermes takasagoensis 371 
Macrotermes barneyi 374 
Reticulitermes flaviceps 470 
Reticulitermes flavipes 470 
Coptotermes gestroi 470 
Coptotermes formosanus 470 
Nasutitermes takasagoensis 465 
Macrotermes barneyi 468 





480 ba 

495 
495 
495 
495 
490 
493 


Reticulitermes flaviceps 
Reticulitermes flavipes 
Coptotermes gestroi 
Coptotermes formosanus 
Nasutitermes takasagoensis 
Macrotermes barneyi 





Fig. 2 Multiple alignments of amino acid sequences of the B-glucosidase genes of 
Reticulitermes flaviceps ( GZRfBG1 ) and other termite species 
Species source and GenBank accession numbers of B-glucosidase genes; Reticulitermes flaviceps ( GZRfBG1 ), AMY26906. 1; Reticulitermes flavipes , 
ADK12988. 1; Coptotermes gestroi, AGS32242. 1; Coptotermes formosanus, ADB23476. 1; Nasutitermes takasagoensis, BAI50023. 1; Macrotermes 
barneyi, AFD33364. 1. The black, dark gray and light gray denote the amino acid sequence identity of 60% , 80% and 100% , respectively. 


the tGZRfBG1 inclusion body, as well as soluble tGZRfBG1 and pET32a-CK in 100 mmol/L SAB (pH 
protein in cell lysed supernatant, and purified soluble 5.6, 30% glycerol) were 510 and 1 080 pg/mL, 
tGZRfBG1. During ultrafiltration of tGZRfBG1, respectively (Fig. 3). 

precipitate was observed at a high protein 3.3 Enzymatic activity of the recombinant PB- 
concentration. The final concentrations of purified glucosidase 
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Fig. 3 SDS-PAGE (A, C) and Western blot (B) analysis of the recombinant GZRfBG1 
M1, M2, M3; Protein molecular weight marker; 1; Total protein of pET32a-CK expressed in E. coli Rosetta 2 (DE3) ; 2; Total protein of tGZRfbg1- 
pET32a plasmid expressed in E. coli BL21; 3; Total protein of tGZRfBG1 expressed in Æ. coli Rosetta 2 (DE3) ; 4: Supernatant of crude cell lysate of 
tGZRfBG1; 5; Purified soluble tGZRfBG1; 6; Total protein of pET32a-CK expressed in Æ. coli Rosetta 2 ( DE3 ); 7: Total protein of tGZRfBG1 
expressed in E. coli Rosetta 2 (DE3) ; 8; Supernatant of crude cell lysate of tGZRfBG1 ; 9; Purified soluble tGZRfBG1; 10; Purified soluble pET32a- 


CK; 11; Total proteins of pET32a-CK expressed in E. coli Rosetta 2 (DE3 ) . 


The pET32a-CK solution was used as a negative 
control. The cellulolytic activity of purified tGZRfBG1 
was observed in the CMC-containing agarose plate. 
The halo diameters formed by 2.04 (4 uL) and 4. 
08 ug (8 uL) purified tGZRfBGI were 5 and 9 mm, 
respectively. No halo was formed by equivalent 
amounts of the pET32a-CK solution (Fig. 4). The 
specific activity of purified tGZRf{BG1 to cellobiose , 
CMC and salicine, measured by reducing sugar 
production and expressed as glucose equivalent, was 
0.609 +1.27E-01, 0.087 +3.41E-02 and 0.009 + 
1. 74E-03 U/mg, respectively (Fig. 5). Errors 
were given as standard deviation of three independent 
experiments. 


4 DISCUSSION 


Termites are the key cellulose consumers and 
digesters in natural ecosystem, and they have the 
outstanding ability to degrade lignocelluloses using 
their unique enzyme complexes (Sun and Scharf, 
2010; Sun, 2014 ). 
representative wood-feeding termite 


China. A full-length cDNA 


R. flaviceps is one of a 
species in 
sequence of B- 





Fig. 4 Enzymatic activity of purified tGZRfBG1 
on an agarose plate containing 0.2% CMC 
Equal amounts ( pg) of individual tGZRfBGI and pET32a-CK 
preparations were spotted on the plate. CK indicates the negative control 
pET32a-CK. The plate was incubated at 37°C for 18 h and the halos, 
indicating activity, were shown after staining with 0.1% Congo red and 


de-stained with 1 mol/L NaCl. 


glucosidase designated as GZRfbg1 was cloned from 
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0.00 
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Enzymatic category 
Fig. 5 Enzymatic activity of purified tGZRfBG1 
toward different substrates 
The specific activity of purified tGZRfBG1 toward cellobiose (CB) , 
CMC (EG) and salicine (BG) was measured by reducing sugar 


production and expressed as glucose equivalent. Each value is the mean 
+ SD of three independent experiments. 


R. flaviceps and heterologously expressed in E. coli 
Rosetta 2 ( DE3 ) . 


enzymes adopted by termites may help biomass 


Understanding the cellulolytic 


utilization projects. 

In 2008, Xue et al. (2008) firstly purified a 
protein from the salivary gland of R. flaviceps with B- 
glucosidase activity. Its nucleotide or amino acid 
sequence remains unknown, but its predicated 
molecular mass is 93.6 kD which is much larger than 
that of GZRfBG1 (56.45 kD). Multiple endogenous 
B-glucosidases with different properties and functions 
had been found in C. formosanus and Reticulitermes 
speratus Kolb (Shimada and Maekawa, 2014; Feng 
et al., 2015), suggesting that different endogenous 
B-glucosidases maybe also encoded by R. flaviceps. 
Based on the amino acid sequence, GZRf{BG1 can be 
categorized as a GHF1 enzyme. It has a digestive 
function similar to other endogenous termite B- 
glucosidases, including RfBGluc-1 from R. flavipes 
(Scharf et al., 2010; Zhang et al., 2012b) , MbEG1 
from Macrotermes barneyi ( Wu et al., 2012), 
GINkBG from N. koshunensis ( Uchima et al., 
2011), CfBGI from C. formosanus ( Zhang et al., 
2010; Zhang et al., 2012a). Purified tGZRfBG1 
exhibited increased specific activity toward salicine, 
CMC and cellobiose, and also formed hydrolyzation 
CMC is normally 
used as the substrate for EG activity measurement 


circles on CMC-agarose plates. 


and salicine is typically used as the substrate for 
measurement of BG activity ( Eveleigh et al., 
2009). The RfBGluc-1 also showed relatively higher 
activity to cellobiose than to other substrates. And 
the activity of B-glucosidase protein of 93. 6 kD 
purified from R. flaviceps exhibited similar optimum 


temperature, pH and substrate preference to 
RfBGluc-1 ( Xue et al., 2008; Scharf et al., 2010). 
Multiple revealed that the 
deduced amino acid sequence of GZRfBG1 shares 
95% identity with RfBGluc-1. These 
suggested that the salivary gland B-glucosidase from 
sibling species R. flaviceps and R. flavipes are very 
conservative in 


alignment analysis 


results 


sequence and hydrolysis 
characteristics. 

The low B-glycosidase activity of tGZRf{BG1 
might because that the insect originated cellulases 
expressed in E. coli could not undergo the complete 
processing or  posttranslational modification that 
occurs in eukaryotic systems. Bioinformatics analysis 


suggested that GZRfBG1 has 


modification allowing secretion into tissues. The Æ. 


post-translational 


coli expression system used here was convenient and 
provided comparability with other research ( Ni 
et al., 2007; Wu et al., 2012). But mass 
production of native endogenous cellulases from 
termites in E. coli can be difficult (Ni et al. , 2005; 
Hirayama et al., 2010; Ni et al., 2014). In this 
study, tGZRf{BG1 was not expressed in E. coli BL21 
but successfully expressed in E. coli Rosetta 2 DE3 
which is specialized for eukaryotic minor codon 
preference. In addition, the IPTG concentration and 
expression temperature were optimized to increase 
the solubility of tGZRfBG1 ( data not shown ). 
of the protein tGZRfBG1 was 


expressed as an inclusion body and precipitated 


However, most 


during ultrafiltration centrifugation. So, the 
concentration of purified tGZRfBG1 was limited and 
this could be another reason for the very low 
cellulase activity detected in purified tGZRfBG1. 
Differences in the expression system or in the 
expression condition can have significant effects on 
GZR{BG1 


endogenously produced termite B-glucosidase in 


expression and purification. Many 


eukaryotic expression systems ( P. pastoris, 
Aspergillus oryzae, and baculovirus ) were more 
thermostable or active than those expressed in E. coli 
(Ni et al. , 2007; Scharf et al., 2010; Uchima et al. 
, 2011; Wu et al., 2012). As the next step, 
additional properties of GZRfBG1 will be studied 
using a eukaryotic expression system. 

This study characterized a B-glucosidase from 
R. flaviceps, which can provide a theoretical basis 
for the development and phylogenetic study of termite 


cellulases in the future. 
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摘要 :【 目的 】 黄 肢 散 和 白蚁 Reticulitermes flaviceps 是 我 国 代表 性 的 食 食 木 白蚁 之 一 ,具有 高 效 的 纤维 素 
酶 系统 。 本 研究 旨 在 克隆 黄 肢 散 白 蚁 中 纤维 素 酶 基因 并 进行 其 特性 分 析 , 为 白蚁 纤维 素 酶 资源 的 
开发 提供 理论 基础 。【 方 法 】 采 用 反 转 录 PCR( 简 并 PCR 以 及 cDNA 末端 快速 扩 增 实验 
(RACE) 法 克隆 黄 肢 散 和 白蚁 B- - 糖 普 酶 基因 序列 并 进行 生物 信息 学 分 析 ; 利 用 大 肠 杆 菌 Escherichia 
coli 原核 表达 该 基因 ,表达 出 的 蛋白 质 经 纯化 后 进行 纤维 素 酶 活力 测定 。【 结果 】 从 黄 肢 散 白 蚁 中 
克隆 到 一 个 B- 糖 普 酶 基因 CZRfbgl (GenBank 登录 号 ; KU886065 ) 的 全 长 cDNA 序列 ,全 长 1 691 
bp, 开放 阅读 框 1 488 bp, 编 码 448 个 氨基 酸 残 基 ,预测 分 子 量 为 56.45 kD ,具有 典型 的 糖 基 水 解 酶 
家 族 1 (CHF1 ) 的 蛋白 功能 结构 域 。 氨 基 酸 序列 多 重 比 对 显示 , CZRfBGC1 与 北美 散 和 白蚁 
Reticulitermes flavipes 的 唾液 腺 B-glucosidase 氨基 酸 序列 一 致 性 为 95%。GZRIBG1 在 大 肠 杆 菌 
Rosetta 2 (DE3) 中 进行 体外 表达 ,经 纯化 的 GZRfBG1 纤维 素 酶 活力 测定 结果 显示 ,其 对 水 杨 昔 、 羧 
甲 基 纤 维 素 钠 (CMC ) 以 及 纤维 二 糖 的 水 解 活性 依次 上 升 , 能 够 在 CMC- 琼 脂 平板 中 产生 清晰 的 水 
FA. 【结论 ] 结果 表明 ,来 自 黄 肢 散 和 白蚁 与 北美 散 和 白蚁 两 个 近 缘 和 白蚁 的 唾液 腺 B- 糖 葫 酶 在 氨基 酸 
序列 及 其 水 解 特征 方面 亦 具 有 高 度 保守 性 。 该 B- 糖 普 酶 基因 的 研究 丰富 了 和 白蚁 纤维 素 酶 开发 及 
系统 发 生 研究 的 理论 基础 。 
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